Black holes in the TeVeS theory of gravity and their thermodynamics 
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TeVeS, a relativistic theory of gravity, was designed to provide a basis for the modified Newtonian 
dynamics. Since TeVeS differs from general relativity (e.g., it has two metrics, an Einstein metric 
and a physical metric), black hole solutions of it would be valuable for a number of endeavors 
ranging from astrophysical modeling to investigations into the interrelation between gravity and 
thermodynamics. Giannios has recently found a TeVeS analogue of the Schwarzschild black hole 
solution. We proceed further with the program by analytically solving the TeVeS equations for a 
static spherically symmetric and asymptotically flat system of electromagnetic and gravity fields. 
We show that one solution is provided by the Reissner-Nordstrom metric as physical metric, the 
TeVeS vector field pointing in the time direction, and a TeVeS scalar field positive everywhere 
(the last feature protects from superluminal propagation of disturbances in the fields). We work out 
black hole thermodynamics in TeVeS using the physical metric; black hole entropy, temperature and 
electric potential turn out to be identical to those in general relativity. We find it inconsistent to 
base thermodynamics on the Einstein metric. In light of this we reconsider the Dubovsky-Sibiryakov 
scenario for violating the second law of thermodynamics in theories with Lorentz symmetry violation. 

PACS numbers: 04.70.-s, 97.60.Lf, 04.70.Bw, 04.70.Dy 
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I. INTRODUCTION 

There are significant discrepancies between the visi- 
ble masses of galaxies and clusters of galaxies and their 
masses as inferred from Newtonian dynamics. In par- 
ticular, the accelerations of stars and gas in the out- 
skirts of galaxies or those of galaxies in clusters are much 
too large, and the disk rotation curves of spiral galaxies, 
which are naively expected to drop as r -1 / 2 away from 
galaxy centers, tend to remain flat to the last optically 
or radio measured point. These discrepancies are also 
manifest in the gravitational lensing by galaxies and clus- 
ters. It is commonly assumed that these problems stem 
from the existence in the said systems of large amounts of 
"dark matter" . For instance, galaxies are assumed to be 
enshrouded in roundish halos of dark matter that domi- 
nate the gravitational fields far from the galaxy centers. 

But the putative dark matter has yet to be identified 
or detected directly. Furthermore dark matter models of 
galaxies require much fine tuning of the dark halo param- 
eters to fit the data, and there are some sharp problems 
outstanding such as the observationally inferred absence 
of cusps in the dark matter density at galaxy centers, 
cusps which are predicted by dark matter cosmogony. 
Thus many have wondered if dark matter is the whole 
story. An alternative, if less orthodox, approach is for- 
malized in the modified Newtonian dynamics (MOND) 



paradigm which proposes that Newtonian gravity 
progressively fails as accelerations drop below a charac- 
teristic scale ao ~ 10 _10 m/s 2 which is typical of galaxy 
outskirts. MOND assumes that for accelerations of order 
ao or well below it, the Newtonian relation a = — V$Ar 
is replaced by 



fi(\a\/a )a = -V$ 



N, 



(1) 
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where the function fl(x) smoothly interpolates between 
fl(x) = x at i < 1 and the Newtonian expectation 
fi(x) = 1 at x 3> 1. This relation with suitable standard 
choice of fi(x) in the intermediate range has proved suc- 
cessful not only in in rationalizing the asymptotical flat- 
ness of galaxy rotation curves where acceleration scales 
are much below ao, but also in explaining detailed shapes 
of rotation curves in the inner parts in terms of the di- 
rectly seen mass, and in giving a precise account of the 
observed Tully-Fisher law which correlates luminosity of 
a disk galaxy with its asymptotic rotational velocity [|J . 

The pristine MOND paradigm does not fulfill the usual 
conservation laws, does not make it clear if the departure 
from Newtonian physics is in the gravity or in the iner- 
tia side of the equation F = ma, and does not teach us 
how to handle gravitational lensing or cosmology in the 
weak acceleration regimes. All these things are done by 
TeVeS [3], a covariant field theory of gravity which has 
MOND as its low velocity, weak accelerations limit, while 
its nonrelativistic strong acceleration limit is Newtonian 
and its relativistic limit is general relativity (GR). TeVeS 
sports two metrics, the "physical" metric on which all 
matter fields propagate, and the Einstein metric which 
interacts with the additional fields in the theory: a time- 
like dynamical vector field, u, a scalar field, (j>, and a 
nondynamical auxiliary scalar field a. The theory also 
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involves a free function !F, a length scale t, and two pos- 
itive dimensionless constants k and K. 

TeVeS is an attempt to recast MOND into a full phys- 
ical theory in which the latter's novel behavior is due to 
the gravitational field. Some checks of its consistency 
and comparisons with hard facts have been made. Thus 
Bekenstein showed that TeVeS's weak acceleration limit 
reproduces MOND, and that it also has a Newtonian 
limit, and calculated its parametrized post-Newtonian 
coefficients (3 and 7,which agree with the results of so- 
lar system tests Jl 0, Q . Skordis et al. [1] and Dodel- 
son and Liguori 7] studied the evolution of homogenous 
and isotropic model universes in TeVeS, and showed that 
it may reproduce key features of the power spectra of 
the cosmic microwave background and the galaxy dis- 
tribution. TeVeS has also been tested against a mul- 
titude of data on gravitational lensing (for some refer- 
ences see Ref. 0) . All the above refer principally to situ- 
ations where the gravitational potential is small on scale 
c 2 . Since neutron stars and black holes exist in nature 
one must also understand strong gravity systems in the 
TeVeS framework. 

A beginning in the investigation of the stronggravity 
regime of TeVeS has been made by Giannios [5i]. For 
vacuum spherically symmetric and static situations he 
showed that under a simplifying limit (which we shall de- 
tail below), the Schwarzschild metric qua physical metric 
and a particular scalar field distribution together consti- 
tute a black hole solution of TeVeSs. This motivates us 
to look in this paper at more complicated cases, such as 
that of the charged nonrotating black hole in TeVeS. We 
find that the Reissner-Nordstrom (RN) metric as phys- 
ical metric and the usual electric field together with a 
special configuration of TeVeS's scalar field constitute a 
black hole solution in TeVeS. Using this solution we in- 
vestigate the thermodynamics of spherical black holes in 
TeVeS. 

In Sec. fll] we recapitulate the fundamentals of TeVeS, 
while in in Sec. Mil we describe Giannios' results for the 
nonrotating vacuum black hole. In Sec. IIVI we go on to 
solve the TeVeS equations for the case of a charged non- 
rotating black hole, obtaining a physical metric which 
coincides with the RN metric of GR. Sec. [V] presents a 
resolution of the problem pointed out by Giannios: the 
uncharged black hole solution he found seems to per- 
mit superluminal propagation near the black hole hori- 
zon. Next in Sec. I VII we examine how the familiar con- 
cepts of black hole thermodynamics apply to our black 
hole solutions, and check their consistency using several 
prescriptions. We calculate the relevant thermodynamic 
quantities using the physical metric, and show that the 
Einstein metric is inappropriate for discussing thermody- 
namics. In this light we discuss anew the potential ther- 
modynamic inconsistency described by Dubovsky and 
Sibiryakov for theories with broken Lorentz symmetry Q . 



II. THE TeVeS EQUATIONS 

The acronym TeVeS refers to the Tensor- Vector-Scalar 
content of the theory. The tensor part pertains to the two 
metrics, t^ t „, dubbed the Einstein metric, on which the 
vector and the scalar fields propagate, and the physical 
metric g^, on which matter, electromagnetic fields, etc. 
propagate. The physical metric is obtained from the Ein- 
stein metric through the following relation: 



9a(3 



-2<t 



g a p - 2u a U/3 sinh(20). 



(2) 



Thus one passes from the space of g a p to that of g a p by 
stretching spacetime along the vector u by a factor e 2 ^, 
and shrinking it by the same factor orthogonally to that 
vector. This prescription retains MOND phenomenology, 
while augmenting the gravitational lensing by clusters 
and galaxies to fit observations. 

The dynamics of the metrics and the fields are deriv- 
able from an action principle. The action in TeVeS is the 
sum of four terms. The first two are the familiar Hilbert- 
Einstein action and the matter action for field variables 
collectively denoted /: 



>S<> 



"5m — 



f 



16ttG 



g a0 R a pV=9<Px, 



g d x. 



(3) 
(4) 



Next comes the vector field's action (K is a dimensionless 
positive coupling constant) 



S v 



K 



32ttG 



(s Q V«[a,/*]U|j8,./]) 



2A 



— (g^u^ + l) V^gd^x, (5) 



which includes a constraint that forces the vector field to 
be timelike (and unit normalized) ; A is the corresponding 
Lagrange multiplier. The presence of a nonzero u a estab- 
lishes a preferred Lorentz frame, thus breaking Lorentz 
symmetry. Finally, we have the scalar's action (k is a 
dimensionless positive parameter while t is a constant 
with the dimensions of length, and T a dimensionless 
free function) 



S, = 



f 



2k 2 e 2 G 



F(kfh a(3 <t> ta (f>, ) 



(6) 



Above h af3 = g af3 - u a u with u a = g^up. The scalar 



action is here written differently than in Ref. [H; we have 
eliminated the nondynamical field a and redefined the 
function T . The new form makes it easier to understand 
the strong acceleration limit of the theory, which is espe- 
cially relevant to the present work. 

Variation of the total action with respect to g a ° yields 
the TeVeS Einstein equations for g a p; 



G, 



aj3 



8ttG T, 



'a/3 



(1 



7 Q( 3 

(7) 
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The sources here are the usual matter energy-momentum 
tensor T a p, the variational derivative of S m with respect 
to g a ° , as well as the energy-momentum tensors for the 
scalar and vector fields: 



JL 

kG 



2k 2 e 2 G ' 



(8) 



7a/3 



with 



\u a Uf3 



fx(x) = T (x). 



(9) 



(10) 



Each choice of T defines a separate TeVeS theory, and 
n(x) is similar in nature to the function jl in MOND. 
In particular, fi(x) ~ 1 corresponds to high acceleration, 
i.e., to the Newtonian limit. 

The equations of motion for the vector and scalar fields 
are, respectively, 



kG [g a P + (1 + e'^) u a uP] f af} , (11) 



:;/3] 



Sir 



+\u a + — nu^, pS T<f>, 



= 8nG (1 - e"**) g^vPf^ . 



(12) 



Additionally, there is the normalization condition on the 
vector field: 



-1. 



(13) 



The lagrange multiplier A can be calculated from the vec- 
tor equation. 



III. NEUTRAL SPHERICAL BLACK HOLES 

In his work on black holes, Giannios [5[ worked in the 
limit fi — ► 1 which also entails T{x) = x. Since we shall 
later work in the same limit, we shall here justify it in 
more detail than he did. Near the horizon of a black 
hole of mass m, the Newtonian acceleration amounts to 
lO 23 (Af0/m)ao- Thus even for the most massive black 
holes suspected (10 10 m Q ) the accelerations are strong 
on scale clq out to at least a million times the gravita- 
tional radius, i.e. well into the asymptotically flat re- 
gion which determines the metric properties. This means 
MOND effects are suppressed while the full complexity 
of the TeVeS equations is still evident. In the said limit, 
and under the assumption that the vector field points in 
the time direction (which has support in the more gen- 
eral context of static solutions [3]), Giannios obtained 
an exact spherically symmetric analytical solution to the 
TeVeS equations, for metric, scalar and vector fields. 

The Einstein metric is taken in isotropic coordinates, 



.0 _ 



t,x L 



9 and 



„3 _ 



ip, 



ds 2 



gc^dx 13 = -e v dt 2 + e c {dr 2 + r 2 dVt 2 ), (14) 



where henceforth dfl 2 = d9 2 + sin 2 9 d6 2 . Since v, £ are 
functions of r only, and the vector field points in the 
time direction, its r dependence is fully determined by 
the normalization condition Eq. (|13[) and the requirement 
that u a be future pointing: 



i a = {e-^ 2 , 0,0,0). 



(15) 



Then the relation between the physical and fields metrics 
reduces to 



gtt = e^gu, 

Qii ~ C Qii 



(16) 
(17) 



Giannios first solved the TeVeS equations assuming 
that K — 0, thus decoupling the vector field from the 
theory, and then performed a transformation involving 
K, which recovered the more general solution. For K = 
the Einstein's tt 1 rr and 99 equations are, respectively, 



r 4 



47r(0' 



l\2 



C + v' , (C) 



'\2 



r '42 

v' + C (y'Y + X" + tv" 

2r + \ 
and the scalar equation takes the form 



k 

4ir((j)') 2 
k ' 
4tt(0') 2 



2r 



= 0. 



(18) 
(19) 
(20) 

(21) 



Since there are only three unknown functions, i / (f),C( r ) 
and 4>[r) , one of the four equations is obviously superflu- 
ous. 

Combining the rr and 99 Einstein equations gives the 
simple differential equation 

2(v + C)" + 6{ " + CY +((^ + C) , ) 2 = 0- (22) 



This has the solution 



2 In 



(23) 



where the additive integration constant has been set to 
zero in order to have an asymptotically flat spacetime, 
namely, v, C, — ► when r — > oo. 

The second integration constant, r c , can be evaluated 
by expanding v + £ above in 1/r and comparing with the 
1 jr expansions (with K = 0) of the metric coefficients of 
the exterior solution for a spherical mass |3, 5j], 



r 1 r 

r lb 7r 



Gm, 



Here m s is a mass scale 0] defined by the expansion 

kGm s 



(r) = 0c - 



4nr 



(24) 
•(25) 

L 

(26) 
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for the solution of Eq. (|2Tj) . For a ball of nonrelativistic 
fluid, m s is very close to the Newtonian mass, and r g is a 
scale of length that can be linked to the object's mass [|[. 
However, the relation between Gm s and r g depends on 
the system under consideration, and is different for stars 
and black holes. At any rate, for K = 0, r c is found to 
be 



(27) 




Making the educated guess that 
4r? 



C' = 



r(r 2 — r 2 ) r 2 



Giannios determines v to be 

v = —Z- In 

2r r 



(28) 



(29) 



The correctness of Eqs. (|28|) and (J29J) can be checked by 
substituting them into the sum of Eq. (fl~8|) and Eq. (fT9|) . 
or the difference of Eq. (fT5|) and Eq. (|20p ; both combina- 
tions are independent of the equation pair already used. 
The determination of the Einstein metric for K = is 
completed by the trivial integration of Eq. (f2"8")l . Finally, 
the scalar field is found now by integrating Eq. (f2"Tj) and 
fixing the two integration constants just as in Eq. (|26p . 



kGm s 

87TTV 



In 



(30) 



Going on to the more general case K ^ 0, Giannios 
finds that just by replacing Eq. (|2T)l by 



(31) 



in the above solutions for i>, C, and 4> will produce an exact 
solution of the TeVeS equations for K ^ [equations 
which are the Q = case of Eqs. (|4"9] ) - (|52]) below]. 
The physical metric now follows from Eqs. (fT6 |) -(fl~7 |l : 




9tt 



(r 2 



(32) 

a 

(33) 

( _ ( . In order for this result to represent 
a black hole, the candidate event horizon r = r c must 
have bounded surface area, and must not be a singu- 
lar surface. The surface area is proportional to g r r{r c ), 
which has a factor (r — r c ) 2 ~ a ; for this to be bounded 
requires a < 2. The Ricci scalar of the above metric is 



with a 



kGm s 



R 



2(a 2 - 4)r 2 r 4 (r - r c ) a ~ 



r c ) a 



+4 



(34) 



We notice that R will blow up as r — > r c unless a = 2 or 
a > 4. Thus, only the value a = 2 describes a black hole. 
The definition of a then gives another relation between 
ry and r„ 



fcGm, 



" 4 8vr ' 
and the physical metric takes the final form 



9tt 



(35) 

(36) 
(37) 



which we recognize as the Schwarzschild metric in 
isotropic coordinates. 

Unlike GR's Schwarzschild black hole, the TeVeS neu- 
tral spherical black hole is "dressed" with a scalar field <fi, 
a solution of Eq. (f3T)|) . This field does not induce a singu- 
larity at the horizon because of the particular structure 
of the TeVeS equations. However, the logarithmic diver- 
gence of <p at the horizon was a cause of concern to Gian- 
nios. It was earlier shown Q that absence of super lumi- 
nal propagation of the various TeVeS fields is guaranteed 
only when <fi > 0. But here cf> diverges logarithmically at 
r = r c , and becomes already negative sufficiently close to 
r c even if <f> c > 0. We will show in the next section how 
this apparent problem is solved. 



IV. CHARGED SPHERICAL BLACK HOLES 

The next natural step is to look for an electrovacuum 
static spherically symmetric solution to the TeVeS equa- 
tions, the analog of the RN solution of GR. We again 
take fi w 1. Again we assume that the vector field points 
in the time direction, and that both the physical and the 
Einstein metrics are spherically symmetric. These are es- 
sential simplifying assumptions which enable us to find a 
specific solution to the TeVeS field equations. Other solu- 
tions may exist for which the vector field is endowed with 
a radial component. However, to judge from the neutral 
case, as analyzed by Giannios [5j, the PPN parameter j3 
of such a solution with very low charge would be in con- 
tradiction with recent observations [9j in the solar sys- 
tem. It would be odd if the PPN structure of a black 
hole's far field were that different from the sun's. By 
contrast, still in the neutral CclSG, cL TeVes solution with 
the vector field pointing in the cosmological time direc- 
tion yields PPN parameters identical to those of GR 0] . 

We continue to work in isotropic coordinates, for which 
the transition between physical and Einstein metrics is 
simplest: as seen earlier, in view of the the normalization 
condition (fl"3")) . the transformation ^ is equivalent to 



9a(3 



e 2 +9tt ■ 



r,U,(j 



(38) 



5 



The Einstein metric again takes the form (| 14[) . and the 
physical metric will have similar form, namely 

d~s 2 = g ap dx a dx p = -e 9 dt 2 + e c (dr 2 + r 2 d£l 2 ), (39) 

with the following relation among i>, £, v and £: 



C(r) =C(r)+2^(r), 
i/(r) = P(r) — 2<f){r). 



(40) 
(41) 



Now, the cosmological value of should be nonzero 
in our evolving universe: cj>(r — > oo) = </> c . Thus the 
requirement that the Einstein metric be asymptotically 
Minkowski (both £ and i/ vanish as r — > oo), needed to 
maintain consistency with previous work [3j, introduces 
a factor ±2</> c in the physical metric coefficients, 



C(r 
v(r 



oo ) 
oo) 



-2<p c 
2(f> c . 



(42) 
(43) 



This is equivalent to a rescaling of the coordinates which 
depends on cosmological epoch, and will have to be taken 
into account when considering physical quantities in the 
framework of TeVeS. 

The energy-momentum tensor no longer vanishes; it is 
given by 

f c0 = ^ (F ap Fp " - ±~g a pF pa F^ , (44) 

with F a f3, the electromagnetic field tensor (not its dual), 
obtained by solving Maxwell's equations in vacuum writ- 
ten wholly with the metric g a p, namely, 



a/3 _ 



(-sT 1/! »f[(-j) l W^ 



= 0. (45) 



In the isotropic metric Eq. (|39j) . and with the assump- 
tion of spherical symmetry and absence of magnetic 
monopoles, the only nonvanishing component of the elec- 
tromagnetic field tensor is 



(46) 



The constant of integration Q will be shown in Sec. IIVI 
to coincide with the physical electric charge of the black 
hole. 

Since we assumed the vector field to point in the time 
direction, then as in the vacuum case, the normalization 
condition (|13|) determines its functional dependence: 



u a = {e~ v ' 2 , 0,0,0). 



(47) 



It follows that the spatial components of the vector equa- 
tion (fl"2"|) are identically satisfied, while its temporal com- 
ponent serves to determine the Lagrange multiplier A to 
be substituted in the Einstein equations: 



A : 



K {ris'C + Irv" + 4i/) GQ 2 e 2 *(e 4 * - 1) 



4 p 2C 



(48) 



We now turn to the Einstein equations (O, and the 
scalar equation (fTTj) . Upon substitution of the Lagrange 
multiplier |48|) and the electromagnetic field tensor (|46p , 
the tt 1 rr and 99 equations become, respectively, 

2C (C) 2 // K (%v' + 2rv'C + r{v') 2 + 4rv") 



47r(</>') 



IV2 



GQ 2 



k r 4 

C+jS (CO 2 CV Kju'f 

r 4 2 8 

4^(f) 2 e-^GQ 2 



(v' + C) + 2C" + 2i/") K{v'f 



(49) 



(50) 



2r 



4ar(<p 



>\2 



-C+24, 



GQ 2 



The scalar equation is 

jj, (r W + O + 4) ft er^ kGQ 2 



2r 



47rr 4 



(51) 



(52) 



These are four equations for three unknowns C( r )> v i r ) 
and <p(r), so one of the equations is actually redundant. 
We shall use two combinations of the three Einstein equa- 
tions plus the scalar equation. 

By adding the rr and 99 equations, we again obtain 
as in the vacuum case Eq. (|22p. This time we write the 
solution 



C + v = 2 In 



(53) 



Here one integration constant has been set so as to have 
an asymptotically flat spacetime, namely, i/, £ — *■ when 
r — > 00. The other constant, r^, will be set by the bound- 
ary conditions on the horizon. 

The remaining equations for v, £ an d </> are not im- 
mediately solvable. To make progress we shall assume 
that the physical metric g a p is of RN form, solve for the 
scalar field in this framework, and check that all TeVeS 
equations are satisfied. This will give us a pair of charged 
black hole solution of TeVeS; existence of other solutions 
is yet to be excluded. 

In Schwarzschild coordinates x° = t, x 1 = R,x 2 = 9 
and x 3 = <p, the RN metric may be written as 



ds 2 = - (1 - R+/R) (1 - R-/R) dt 2 
dR 2 



(1 - R+/R) (1 - R-/R) 



R 2 dil 2 (54) 



where i?+ and i?_ are the coordinates of the outer and in- 
ner horizons, respectively. We may transform the metric 
to isotropic form by going over to a new radial coordinate 
r defined implicitly by 

R(r) = r+{R + - RJf/lQr + [R+ + RJ)/2. (55) 
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This gives 



ds 2 



(4r - (R+ - R-)) 2 (4r + (R+ - R_)) 2 
' (16r 2 + 8r(R+ + R-) + {R+ - R-) 2 ) 2 
(16r 2 + 8r{R + + R_) + (R + - i?_ ^ 2 



dt 2 



256 r 4 



-(<ir 2 +r^ft 2 ). 



We recall that in GR i?+ and i?_ satisfy the relations 
i?+ + i?_ = 2Gniti and R+R- = Gno 2 , where we write 
the gravitational constant as Gn to distinguish it from 
plain G, the coupling constant in TeVeS. We shall here 
assume that the physical metric of TeVeS has the above 
form while leaving the parameters i? + and i?_ to be de- 
termined later. However, the proposed metric is asymp- 
totically Minkowskian, while as previously mentioned, we 
require rather that the Einstein metric be asymptotically 
Minkowskian. This means that in the generic physical 
metric Eq. ([55)) we must set 

- M 1 (±r - (R + - R_)) 2 (4r + (R + - R_)) 2 

V \ r ) = 111 777. 9 , n 1 7-) \ i / n r, N9N9 + 2( Pc, 



(16r 2 + 8r(R+ + R_) + (R+ - i?_) 2 ) 2 



(56) 



- (16r 2 + 8r(i? + +i?_) + (i 1 ' + -i 1 '_) 2 ) 2 
C«=ln ^ 20 c . 

(57) 

To simplify these note that by Eqs. (|4"D")) - (l4"l"]) we have 
C + v = C + i>, whereupon in view of Eq. (|53p . 



(r 2 - r 2 ) 2 _ (4r - (R + - i?_)) 2 (4r + (R + - R_)) 2 
r 4 256 r 4 

(58) 

We may thus relate R + and i?_ to the integration con- 
stant r/j appearing in 



(59) 



r h = -(R+-R-) 

4 



Since R{r = rh) = R + , r = rh is the outer black hole 
horizon in isotropic coordinates. In terms of r/j and M = 
(i?+ + i?-)/2 the physical metric coefficients are 



£ _ (r 2 + r 2 , + Mr) 2 
(r 2 +r 2 + Mr) 2 



(60) 
(61) 



It is useful at this point to trade the charge Q for a 
dimensionless positive parameter a defined by 



This replaces the GR relation = GnQ 2 ■ The value 

of a will be determined by the Einstein equations (|49|) - 

(HQD. 



The only indeterminate function remaining now is the 
scalar field. The scalar equation (|52|) can be rewritten in 
terms of the new parameters M, rh and a as 



2r4>' ka 2 (M 2 - 4r 2 > 



2<*,. 



Its general solution is 



fee 2 * 1 



47r(r 2 



Mr) 5 



0. 



(63) 



47T 

(1-C) ln(r-r h ) 



x [(1 + C) ln(r + r h ) (64) 



ln(r 2 



Mr)], 



with </> c and C integration constants, the first already fa- 
miliar. Since we guessed the form of the metric, we need 
to verify that the Einstein equations are satisfied. From 
the requirement that Eq. ([50]) be satisfied, we obtain val- 
ues for a and C: 



C± 



4tt(2 - K)e~ 2 ^ 
k(2 - K ) + 8tt ' 
y / 2k 2 (2 — K) + 8nkK 
{2-K)k 



± 



(65) 
(66) 



Eq. (|4"9")) is then satisfied identically. Since we have al- 
ready used the sum of Eq. ([50]) and Eq. (fBTj) to get the 
solution ([SB")) , we see that all TeVeS equations are satis- 
fied. Thus the RN metric from GR with a suitable choice 
of parameters is the physical metric of TeVeS spherical 
charged black holes. 

We shall soon see that a physically acceptable solution 
can be had only for K < 2. For such solutions the sign 
of the quantity under the square root in Eq. (f6"6"]) is posi- 
tive. The two TeVeS solutions (corresponding to the two 
signs of C) are most clearly presented in terms of the 
coefficients 5± = (k/A-K)a 2 (l + C±)e 2 ^, or 



(2 - K)k ± y/2k 2 (2- K) + 8irkK 
(2-K)k + 8-K 



(67) 



GQ 2 = a 2 R+R- = a 2 (M 2 - 4r/J 



(62) In view of Eq. (|4"T|) we finally obtain the solutions 
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,~2 



cf>(r) 



(r 2_ r 2 )S 



2 ( r2 + r h + Mr f 



(r 2 + r 2 + Mr) 2 
[r + r h ) 

( r 2 +r 2 + Mr)( i ++ { -- 2 )/ 2 



(/) c + 6± ln(r + r h ) + S T ln(r - r h ) - -(£+ + <5_) ln(r 2 + + Mr), 



0,0,0 



(68) 
(69) 
(70) 



RESOLVING THE SUPERLUMINAL 
PARADOX 



We have found two black hole solutions for each value 
of Rh and M. The requirement that superluminal prop- 
agation be excluded selects one of them as physically vi- 
able. As mentioned in Sec. IIII1 in a region where cj> < 
superluminal propagation of the TeVeS fields is not ruled 
out. This acausal behavior would be unacceptable. Now, 
since In (r — rh) in Eq. (|6"9"|) is arbitrarily large and neg- 
ative near enough to the horizon, its coefficient must be 
negative in order that <fr have a chance to be nonnega- 
tive everywhere. It is easy to see that for K > 0, k > 0, 
5 + is always positive. Thus the solution Eqs. (|6"8"j) - (|7Tj|) 
with lower signs is immediately excluded on grounds that 
it permits superluminal propagation. But is the second 
solution viable in this sense? 

Focusing on the solution with upper signs, we must 
now exclude the parameter range K > 2 + 87r/fc; the 
equality here corresponds to unbounded <5_ and </>, while 
the inequality leads to 5— > and superluminal propa- 
gation. The range 2 < K < 2 + 8n/k, although palatable 
in this sense, gives a 2 < 0. We shall show in Sec. lVII that 
this is unphysical. For < K < 2 we have <$_ < while 
a 2 > 0. Thus a physically viable black hole solution of 
TeVeS can exists only for < K < 2 (we continue to 
assume that k > 0). It is the solution with the lower 
indices in Eqs. (|6"8 j) - (|70j) . 

Close enough to the horizon, <j> of this solution is nec- 
essarily positive because of the 5- ln(r — rh) which is 
arbitrarily large. Additionally, the asymptotic value of 
<p is 4> c , the cosmological value of the scalar, which may 
be assumed to be positive Hence the question of 
whether (f>(r) is positive in the intermediate region hinges 
on whether it has a negative minimum outside the hori- 
zon, or not. 

To find out we look at its derivative, 



if 



4>'(r) = 



(M + 2r h )(r + r h ) 2 5_ + (M - 2r h )(r - r h ) 2 6+ 



2(r 2 -r 2 )(r 



+ 4 



Mr) 



The numerator here is quadratic in r and thus has two 
roots. Now in the case K < 2 we have £_ < 0, but 
> 0. Then because M > 2r h [see Eq. ([55]) and the 
following discussion], both roots are real. Furthermore, 



M < 2r h 



5+ - 5. 



■6- 



2r h 



^/2k 2 {2- K) + 8irKk 



(72) 



(2 - K)k 

both roots are at r < 0, so for r > the field <p(r) has 
no minimum and must be everywhere positive. 
Focus now on the case 



M>2r h 



y/2k 2 (2- K) + 8nKk 
{2-K)k 



(73) 



Now (j) does have a minimum outside the horizon. In 
Fig. [T] we plot (j>— 4> c for several values of M/r^. We see 
that unless M/rh is very large, the dip below the axis 
(which grows roughly as hi M/rh) is modest compared 
to unity. Hence, a modest positive <j) c (which is expected 
from cosmological models [3] ) will be enough to keep <j){r) 
positive throughout the black hole exterior, except for 
black holes with exponentially large values of M/rh for 
which a region of negative 4> will occur near the horizon. 

In fact, for = 0, which by Eq. (f59|) means that 
R + = R_, i.e., that the physical metric is extremal RN, 
the two solutions for 6 are identical: 



(2-K)k 



8tt 



ln(l 



(74) 



r 1 {2-K)k 

Thus for K < 2 and r/, = 0, the variable part of <fi is 
negative and can be very large for r <^ M. This is un- 
acceptable as it permits superluminal propagation. We 
may conclude that provided 0<K<2, k>0 and 4> c 
somewhat above zero, the superluminality issue raised 
by Giannios does not arise for the TeVeS charged black 
hole solution with the lower signs in Eqs. (|6"51) - (f7"0"|) . The 
above conclusion does not apply to black holes near the 
extremally charged case. 

What about Giannios' case Q = for which he found 
conditions conducive to superluminal propagation (see 
end of Sec. [m] )? The TeVeS equations (g9])-([52]) are 
smooth with Q, so we may take the limit Q — > of 
their solutions, Eqs. (f6"8"|) - (fT0"]) . In this limit according to 
Eq. (HU), M = 2r h while by Eq. r h = #+/ 4 - Thus 

our metric ((68)) reduces exactly to Giannios' Eqs. (|36|) - 
(|37[) with the obvious identification r c — r^; that is we 
recover the fact that the physical metric is Schwarzschild. 
In the same limit our scalar field solutions (|69|) reduce to 
the pair 



5 T In 



rh 



rh 



(75) 
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FIG. 1: Our solution for <j> — (f> c as a function of r for several 
values of Al/r^, the higher M/r^, the lower the curve. Both 
axes are in arbitrary units. 



whereas Giannios obtained only one scalar solution. We 
notice that the solution with upper sign has (f> > for all 
r > rh provided that we stick to the parameter ranges 
< K < 2, k > for which <5_ < 0. The solution with 
the lower sign has 4> < sufficiently near r = r^; this is 
Giannios' solution, and it is indeed excluded because it 
allows superluminal propagation. 

To sum up, in our study of spherical static black holes 
in TeVeS, we have found a viable charged black hole so- 
lution for the parameter range < K < 2, k > 0. The 
limiting case Q — > of this is a viable neutral black 
hole solution. Since black holes are seen in nature with 
virtual certainty, the above results tell us that only the 
range 0<X<2,fc>0of TeVeS need be considered 
as physical. This range includes the values that have 
been explored in the confrontation of TeVeS with obser- 
vations [1, H, H, • 



VI. BLACK HOLE THERMODYNAMICS 

It has been clear for long that black holes are re- 
ally thermodynamic systems characterized by tempera- 
ture and entropy. Thus a discussion of black hole solu- 
tions in TeVeS would be incomplete without a survey of 
their thermodynamic properties. However, before we can 
talk about thermodynamics for the charged black hole in 
TeVeS, we must first identify the physical values of at- 
tributes of the black hole solution. By physical values we 
mean the quantities than an asymptotically Minkowski 
observer would measure using instruments made of mat- 
ter, measurements which are thus referred to the physical 
metric. These values need not be identical with those of 
quantities naively associated with the attributes. For ex- 
ample, we do not know a priori that the masslike quan- 
tity M and the chargelike quantity Q appearing in our 
solution are indeed the physical mass and charge of the 
black hole. In fact, we shall sec that M is related to 



physical mass in a nontrivial fashion. 

We first note that the G appearing in the TeVeS 
equations is not Newton's constant, but, as shown else- 
where [l(|, is related to it through 



Gn - { 4n(2-K) U " 



(76) 



It will be useful to also write the above relation in terms 
of the constant a defined by Eq. (f6"5| : 



G N = {G/a 2 )e~ 2 ^. 



(77) 



Experimentally Gn > 0; it also seems natural that the 
fundamental coupling constant G be positive; hence we 
must require a 2 > 0. This explains why in Sec. IIVI we 
ruled out the parameter range 2 < K < 2 + 87r/fc. 

Next, recall that if we use the same coordinates for 
the Einstein and physical metric, the transformation ((2|) 
implies that our physical metric is not asymptotically 
Minkowski. Thus, asymptotically, the relation between 
physical distance x and the corresponding spatial length- 
like coordinate (denoted x) must be 

x = e-^x. (78) 

Focus now on M . According to Eq. (|6ip . we may write 
the asymptotical expansion for e u as 



1 



2M 

r 



O 



1 



(79) 



Thus r here is not physical distance f, but is related to 
it through Eq. (|T8[) . Rewriting e v in terms of the latter 
gives 



1 - 



2Me-^ 



(80) 



From the customary linear approximation we see that M 
and physical mass m are related by 



Me" 



GNtn . 



(81) 



Likewise, the physical charge q can be easily identified 
by integrating the flux of the electromagnetic field tensor 
through a spherical shell at spatial infinity: 



1 



q= lim — / F* r eV sin0 d6d(f>. 

r—^oo 47T 



(82) 



Use of in forming the area element guarantees that we 
are calculating a physical flux: according to Eq. (|40| the 
factor e~ 2 ^ c required by Eq. (|75|) is supplied by the e^. 
Substituting F tr from Eq. (|46|) gives 



Thus our charged black hole is characterized by mass m 
and charge Q as measured by physical asymptotic ob- 
servers for which the metric is g a p. 

In investigating the black hole entropy we start with 
the assumption that it is given in terms of the physical 
surface area of the outer horizon A by the usual formula 



S 



A 



BH 



N 



(84) 



It is true that more complicated forms are known, 
but they usually appear in gravity theories with higher 
derivatives; TeVeS is free of these. The proof that 
our choice is correct ultimately rests on the consistency 
checks we present later in this section. 
Obviously 

A = inrle^ = 4tt (2r h + M) 2 e' 2 ^. (85) 
From (|62p we have for the outer horizon 
1 



r h = -^M 2 -GQ 2 /a 2 . 



(86) 



Thus 



A = 4tt (m + v/Af 2 - GQ 2 / 



a 2 ) 2 e' 2 ^ (87) 



We now express A in terms of physical mass, charge 
and Newton's constant using the relations (f8"Tj) , ([77]) : 

A = 4vr (^G N me^ + ^J(G N me^) 2 - G N e 2 ^Q 2 ^j e" 2 ^, 



so that 

Sbh = t^t f G N m + ^(G N m) 2 - G N Q 2 ) * . (89) 

This is identical to the familiar expression for the entropy 
of a RN black hole. To it corresponds the thermodynamic 
temperature Tbh = (9Sbh /Smjg 1 , or 



Tbh — 



h 



y/(G N m) 2 - G N Q 2 



2tt 



(G N m + ^/(G N m) 2 - G N Q- 



(90) 



To check the consistency of our scheme, we now also 
calculate the temperature corresponding to our black 
hole solution using the Euclidean path integral ap- 
proach This approach entails performing a Wick 
rotation of the time coordinate to obtain a Euclidean 
metric. The path integral for the gravitational action 



then becomes the partition function for a canonical en- 
semble. Regularity of the new coordinate system near the 
horizon requires the new time coordinate to be periodic, 
and the period is related to the black hole temperature: 
T = ?i/period. 

We first define I, the radial proper distance from the 
horizon using the physical metric 



all = 



_ (r 2 +r 2 +Mr) , 



dr , 



I = {r-rl/r + M\n{r/r h ))e-' t ' c 



(91) 



Consequently the physical (2-D) line element da 2 for 
fixed 6 and 4> following from metric (|68p becomes 



dd 2 = - 



(r 2 -r 2 ) 



2\2 



3 2 ^dt 2 



(r 2 + r 2 + Mr) 2 ' 
Near the horizon, where r ~ rh, we have 

I w (2 + M/r h )e-^(r - r h ) + O ((r - r,,) 2 ) 



(92) 



(93) 



Substituting this into da 2 and replacing e^ c dt, the global 
physical time interval according to Eq. (|80p . by idr, we 
obtain an expression for the Euclidean metric near the 
horizon 



2 _ ( 2r h e*° 
_ {(2r h + M) 2 



l 2 {dr) 2 +df 



(94) 



For this metric to be regular at I = (r = rh), we 
must regard r as an angular variable with period 



2tt 



(2r h + Mf 



2r h e*" 1 
the corresponding temperature is thus 

Thrue^ 
ir{2r h + M) 2 ' 



Tbh = 



(95) 



(96) 



By means of Eqs. (|86l . (|77| and (|8T|) this can be reduced 
to precisely the form (|90|) . Thus far the thermodynamic 
description based on Eq. ([54^1 is consistent. 

Of course, our black hole solution must exhibit an 
elect ric p otential. By thermodynamics we would expect 
that [if $ M = ~T BH (dS BH /dQ) m - This gives 



BH 



Q 



G N m+^/{G N m) 2 -G N Q 2 



(97) 



which agrees with the potential of the RN black hole 
in GR. To verify this result we shall also calculate the 
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electric potential by a strictly mechanical approach using 
the conservation of energy. We expect the increase in the 
black hole's energy due to the fall into it of a charged 
particle to equal the particle's conserved (kinetic plus 
electric potential) energy. 

We set out from the Lagrangian for a charged particle 
with mass [i and charge e 



di 



, dx a dxf 3 ~ dx a 

C = e-^[-^l- 9aP — — + eA a — 



(98) 



Asymptotically g u — > e 2 "^, so that P t — > \ie^ 
eAte^^"; now since the physical global time is given by 
t = e^ c i, the first term in P t is recognized as minus the 
physical value of the particle's rest plus kinetic energy. 
This justifies the normalization we selected for the La- 
grangian. We can thus identify — eA t e~* c , as the phys- 
ical value of the particle's electric energy. The physical 
electric potential of the black hole is inferred from this 
last energy at the horizon, §bh = —A t (rh)e^^ a . 



where A a is the potential for F a p, The e~^ c normal- 
ization will soon be justified. The Lagrangian does not 
depend on t; therefore, we have the conserved canonical 
momentum 



Pi 



dC 
W 

dr 



-6 ~ dt 
Me 9tt^ + 



eA t e 



(99) 



We calculate At by integrating the electric field |46|) 
from infinity to r/, 



A t (r h ) 



F tr dr 



Qe*(*-0 



dr 



Qe 



2r h + M 



(100) 



r 



Substituting here the expression (|86[) for in terms of 
M and Q, and switching to physical mass using Eq. (f8"Tj) . 
we finally get 



A t (r h ) 



G N m+ y/(G N m)* - G N Q* 



(101) 



But we found the physical black hole electric potential 
to be —At(rh)e~^ a 1 so the present method of calculation 
gives exactly the same result, Eq. (|97|) . as the thermody- 
namic computation. 

The above calculations serve as a consistency check 
of the physical values for mass and charge which we at- 
tributed to the black hole. They also demonstrate the 
physical consistency of a thermodynamical description 
of spherical black holes in TeVeS. In particular, they jus- 
tify our guess (|84|) for the form of the black hole entropy, 
and verify that the first law of black hole thermodynam- 
ics [IH[l3[ holds for the TeVeS spherical black holes. All 
this is accomplished by referring all physics to the phys- 
ical metric. However, there is one issue for which one 
must consider the role of the Einstein metric. 

Recently Dubovsky and Sibiryakov (DS) Q showed 
that in a theory with Lorentz symmetry breaking via 
a time-dependent scalar field, in which there is more 
than one maximal propagation speed, it would be the- 
oretically possible to construct a perpetuum mobile that 
would transfer heat from a colder to a hotter region. This 
would be accomplished, via Hawking radiation, by ex- 
ploiting the different temperatures of the nested horizons 
corresponding to massless fields with different propaga- 
tion speeds. 



DS consider a static spherical black hole, and two 
fields, tpi and ip2, which do not interact with each other 
except through gravity, and propagate at different speeds 
ci and C2- Consequently there exist two distinct horizons, 
one for field tpi that radiates a la Hawking with temper- 
ature T\ , and the second for ip2 radiating at temperature 
T2. It is assumed that C2 > ci; the model then gives 
T% > T±. DS assume the black hole is surrounded by 
two nested shells, shell A which interacts with ipi but is 
transparent to ip2, and shell B which interacts only with 
ip2- Shell A has temperature Ta and shell B is hotter 
at Tb- It is also assumed that Ta > T\ and Tb < T<z. 
DS make the innocuous assumption that heat flows from 
higher to lower temperature, with the heat flow increas- 
ing monotonically with temperature difference, and van- 
ishing only when the two temperatures are equal. Then 
they point out that heat will flow from A to the black 
hole via quanta of tpi and from the black hole to B via 
ip2 particles. It is possible to adjust the shell tempera- 
tures so that the two mentioned flows become equal, in 
which case the black hole is in steady state. Then the 
only overall effect is heat flow from A to B, that is from 
cold to hot. The second law thus appears to be violated. 

TeVeS also breaks local Lorentz symmetry, albeit by a 
different mechanism: it is equipped with a timelike non- 
vanishing future-pointing vector field. Further, TeVeS 
possesses two metrics; this feature implies different prop- 
agation velocities for light and for gravitational waves. 
Does the second law fail in TeVeS within some version 
of the DS scenario? In order to construct the appro- 
priate version, we evidently first need to identify the 
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distinct horizons associated with light and with gravita- 
tional waves. In TeVeS light propagates on the null cone 
of the physical metric g a p, and it is evident immediately 
from Eq. that the horizon for light is at r = r/,. 
By contrast, tensor gravitational waves propagate on the 
null cone of the Einstein metric g a p |3[ . 

With (j> given by the upper sign alternative of Eq. 
the Einstein metric corresponding to the physical metric 



9tt =~e v = — . , , _ 2 ■ 7737^ -fi , (102) 



J+24> = (r-r h ) 2S -(r + r h ) 2S + 

r 4( r 2 _|_ r 2 + Mr )(5_+5 + -2) 



(103) 



For this metric the horizon (for gravitational waves) 
could only be at r = r/,, the same location as the hori- 
zon for light! However, as measured with respect to g a p 
that surface's area diverges: since <5_ < for the phys- 
ical K, k region, g rr blows up at horizon, and since the 
metric is isotropic, this alone causes the area of the sur- 
face r — Th to blow up. The curvature scalars calculated 
with g a fj also diverge at , revealing this location to be 
an essential singularity of the Einstein metric, and not a 
horizon. Thus we are in no position to form a well de- 
fined entropy while working in the geometry perceived by 
gravitational waves. Likewise, we cannot obtain a black 
hole temperature: no entropy, no thermodynamic tem- 
perature. A similar problem arises in trying to calculate 
the temperature by applying the Euclidean prescription 
to the Einstein metric: the gtt does not behave like I 2 . 

The appearance of a singularity of g a p at the same 
surface as the physical metric's horizon does not pose an 
insurmountable problem. It was shown by Zlosnik, Fer- 
reira and Starkman [l4| that TeVeS can be reformulated 
as a vector-tensor theory, with a single metric, the phys- 
ical metric. Probably the formal failure to bring out a 
thermodynamics in the Einstein metric reflects the fact 
that TeVeS is at the bottom a one-metric theory, with 
the Einstein metric being no more than a mathematical 
convenience. 

The absence of a thermodynamic temperature for grav- 
itational waves in the TeVeS black hole background most 
likely means that any Hawking-like emission of these 
waves is not thermal. The attempt, a la DS, to iden- 
tify two distinct black hole temperatures for the same 
black hole, each tied to a different maximal propagation 
velocity, thus fails. However, it has been suggested to us 
that the paradox can still arise as follows. One associates 
with the black hole a graviton effective temperature r in 
lieu of, say, T2 in the DS scenario. This r is defined as the 
temperature that the shell B, which interacts solely with 
gravitons, would have to possess in order to just balance 
the black hole's emission power in gravitons (though of 
course without any pretense of detailed balance). Now 
suppose r exceeds X\, the black hole's photon tempera- 
ture. By suitably adjusting Ta and Tb while observing 
the ordering r > Tb > Ta > 7i, it should be possible to 



annul the overall energy gain of the black hole. Then we 
have a DS scenario where heat flows from the low tem- 
perature Ta to the higher Tb without any other change 
taking place. Of course if r < T\, we have to arrange 
things with t < Tb < Ta < T\ to get flow from Tb to 
the higher Ta- We emphasize that with the introduction 
of the effective temperature r, the propagation veloci- 
ties and Lorentz symmetry violation no longer play the 
dominant role they played in the original DS argument. 

In both of the above setups we seem to have a violation 
of the ordinary second law of thermodynamics. This can 
be avoided only if necessarily t — T l7 when the above 
schemes require Ta — Tb so that no violation is possible. 
Our methods in this paper are not suitable for the study 
of Hawking-like radiation, so we cannot here affirm or 
exclude this possibility. But it seems a fair conjecture 
that, in fact, t = T±. 

Following the DS paper Eling, Foster, et al. [H| sug- 
gested a classical mechanism for violating the second 
law in its generalized form within a gravitational theory 
equipped with a timelike vector field that causes Lorentz 
symmetry violation. This mechanism also relies on two 
maximal propagation speeds, with ca > cb- It is imple- 
mented in a spherically symmetric static situation around 
a black hole. The scenario envisages a particle of type 
A and one of type B. They both fall through the hori- 
zon for particles of type A, and interact with each other 
in a region still outside the horizon for type B particles, 
where the time Killing vector is spacelike with respect to 
the metric sensed by A particles. It is then possible for 
the interaction to make the A particle acquire negative 
energy while the B particle's energy remains positive. 
If now A falls through the B particle's horizon, and B 
escapes passing on its way out through ^4's horizon, the 
black hole's mass has been lowered in the process. Under 
mild assumptions, this means the black hole entropy de- 
creased. But since the particle state may have been pure 
all along, there is no ordinary entropy to compensate, 
and so the generalized second law is violated. 

We point out that in our black hole solution there 
is no intermediate region still outside the second hori- 
zon where the Killing vector already has positive norm. 
Eq. (|68|) shows that in the region r < the Killing vec- 
tor is indeed spacelike with respect to the physical metric. 
However, this region is not outside the null surface [also 
r = rh but in Einstein metric Eq. (|102[) ] that might have 
been construed as the horizon for gravitons. Thus nega- 
tive energy particles can be created, but only inside the 
black hole, and the scenario envisaged in Ref. [l5| cannot 
be enacted. 



VII. CONCLUSIONS AND SUMMARY 

We have here derived a pair of charged spherical static 
black hole solutions of TeVeS that, as far as the phys- 
ical metric is concerned, resemble the RN solution of 
GR. The new features are the TeVeS vector field which 
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points in the time direction, and the spherically symmet- 
ric TeVeS scalar field. We have shown that for a wide 
range of TeVeS parameters, the scalar field of one so- 
lution is positive everywhere as long as the field has a 
modest positive cosmological value. This insures that 
superluminal propagation does not takes place in that 
solution's background. Regarding the TeVeS analogue of 
Schwarzschild's black hole earlier exhibited by Giannios, 
we showed that there are actually two separate solutions 
here too. For one of them the evident positivity of the 
scalar field precludes superluminal propagation. This sin- 
gles out the physical solution. 

An element of guesswork entered in both Giannos' and 
our derivations. He had to guess Eq. f28]): we had to 
guess the RN form of the physical metric. Thus in both 
cases it is not clear if the black hole solutions found are 
the unique ones. Proof of uniqueness in both cases is still 
at large. 

By expressing the parameters of the charged black hole 
in terms of physical attributes measurable by material 
Minkowski observers, we calculated the entropy, temper- 
ature and electric potential characterizing the black hole. 
They turn out to be the same as for GR's RN black hole. 



Black hole entropy and temperature cannot be consis- 
tently defined for the Einstein metric, which would have 
been the correct framework for studying Hawking emis- 
sion of gravitational waves. We consider in this context 
a modified version of the Dubovsky-Sibiryakov [8| sce- 
nario for bringing about a violation of the second law 
of thermodynamics out of Lorentz symmetry breaking. 
This violation of the second law can be forestalled if a 
conjectured equality of the effective graviton radiation 
temperature and the photon Hawking temperature holds. 
The scenario described by Eling, Foster, et al. [l5[ for 
bringing about classical violations of the second law in 
theories with Lorentz symmetry violation cannot be be 
implemented with our TeVeS black hole solutions. 
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